Continuous local model for two-dimensional spin ice 
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We propose a classical model Hamiltonian with a ground state presenting a spin ice structure. 
We analyze the introduction of met ast able excitations on this ground state, showing the emergence 
of pairs of magnetic monopoles. The interaction between monopoles and dipoles in the system 
is studied. As a consequence, we obtain an effective nonlocal interaction between monopoles and 
dipoles from a local classical spin model. 
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The emergence of magnetic monopoles in geometrically 
frustrated magnetic materials (spin ice) [li] is one of the 
most surprising and interesting phenomena recently dis- 
covered, raising a great interest to both experimental [2- 
Q and theoretical jsl-llH research. Spin ice are mag- 
netic materials with ferromagnetic interactions whose 
anisotropy jeopardize the alignment of all the spins, re- 
sulting in (degenerate) ground states with zero magneti- 
zation [12]. These ground states obey the so called ice 
rule, with the spins arranged in such a way that, on each 
vertex of the lattice, two of them point outward and the 
other two point inward the vertex, resulting in zero mag- 
netic flux. Recently, it was predicted that excitations 
on the ground state of the spin ice, corresponding to 
the violation of the ice rule by spin inversion, leads to 
the emergence of quasi-particles with magnetic properties 
that justifies their denomination as magnetic monopoles 

i- 

There are some characteristics of the spin ice sys- 
tems that are essential for the emergence of magnetic 
monopoles in both three-dimensional (3D) real spin ice 
compounds, such as Dy2Ti20j or Ho2Ti20j^ or in artifi- 
cial two-dimensional (2D) spi n ice, as 2D square lattices 
of permalloy nano-islands [13[ - the ice rule (divergence- 
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FIG. 1: Vector plot of the ground state (a) and excited state 
(b) of the model Hamiltonian. 



less of the ground state, i.e., no magnetic charges), pe- 
riodicity (of the pyrochlore and artificial square lattices) 
and anisotropy (the magnetic moments are constrained 
to point along fix axes). Also in magnetic charged spin 
ice [4] (honeycomb spin ice), the last two characteristics 
are present. In both 3D and 2D spin ice, the magnetic 
monopoles are joined by a string, a sequence of inverted 
spins. In 3D systems, the amount of energy necessary 
to separate apart two monopoles is finite, leading to free 
monopoles, whereas in the 2D systems the monopoles 
are confined, i.e., the string energy grows with the string 
length [9]. 

The thermodynamics properties of spin ice structures, 
as well as the interaction between monopoles, has been 
successfully described by discrete Hamiltonians that take 
into account the local nearest neighbor interactions and 
non local (long range) dipolar interaction between the 
spins of the system [H, IqI, 0, [i5|. Also, the magneti- 
zation dynamics in artificial spin ice under the action 
of an external magnetic field was recently studied in 
a discrete mesoscopic model[16]. Here, we propose a 
new classical local model Hamiltonian for a continuous 
spin system inspired in the artificial two-dimensional spin 
ice. The ground state of the model can be mapped 
in the discrete artificial spin ice ground state. We will 
show that, as in the discrete case, the introduction of 
metastable excited states on the system causes the emer- 
gence of quasi-particles that can be identified with the 
magnetic monopoles. We will analyze the stability of the 
excitations, the emergence of magnetic monopoles and 
the interactions (action at a distance) between magnetic 
monopoles or dipoles in the system. 

The model for a two-dimensional square latticed spin 
ice can be constructed from the ground state given by 
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where a is the lattice spacing. Throughout this paper, 
we win set a = tt, without loss of generality. Figure 1(a) 
shows the mapping of the spin field, Eq.([T]). We can see 
that at each vertex of the lattice {x = n-zr, y = m7r, n 
and m integer) the ice rule is obeyed: each vertex has no 
magnetic charge (zero magnetic flux). The ground state 
satisfies the relations V.6'o = and V^6'o = — 26'o. 
By introducing the auxiliary field B = ^{-q^ ^), 



we define our model by the Hamiltonian 
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where n^ A and ji are free parameters of the model. From 
now on, we will assume a>: = 1, implying that H^ X and fi 
are given in units of k,. 

Since all terms of the Hamiltonian are positive-definite, 
the minimum of H occurs only if each quadratic term on 
its integrand vanishes. So, the two first terms on the inte- 
grand on Eq.Q implies that the ground state presents di- 
vergenceless and periodicity, as desired. Due to the pres- 
ence of the product of derivatives ^^ -^ on the third 
term, the Hamiltonian is not invariant under continuous 
rotations in the xy plane, although it is symmetric under 
discrete 7r/2 rotations, a feature that is shared by the ac- 
tual artificial 2D spin ice system. This is the continuous 
version of the anisotropy of the spin ice. The presence 
of this term is crucial to the model, since it fixes the 
amplitude of the ground state, so that solutions with low 
amplitudes do not imply, necessarily, in lower energies for 
the system. The ground state of Eq.(|2]) is infinitely de- 
generate, since any translation or discrete 7r/2 rotations 
on So imply also in a solution with zero energy. This is 
the manifestation of the spin ice residual entropy of the 
ground state in the continuous model. 

Next, we introduce an excitation by adding to the 
ground state a localized function that "fiips" the spins 
located around certain position rg- We propose (,r^{r) = 
^-i{r-re) j(^(j^_f^Y)^ where the Gaussian function indi- 
cates that the excitation is localized, with 7 being the 
excitation width, and /((r — re)^) = A{a{f — feY — 
2){cos{0)i + sin{0)j)^ a function centered in f^- The 
free parameters A, a, and 7 are to be fixed in order 
to minimize H. We will refer to A as the amplitude of 
the excitation, A = implying that the system is in the 
ground state. The parameter gives the direction of the 
excitation, and was introduced in order to show that the 
excitation acts in the opposite direction of the ground 
state, locally fiipping the spin. The functional form of 
the excitation ^^^ we employ here was chosen in order 
to contain the fewest possible free parameters, and it 
was tested against other functional forms with the same 
number of parameters. The fourth order power on / was 
also tested against other even powers of |r — re|, and for 
the values of A and /i used here, showed up to be the 
best choice. Besides, we also found out that there are no 
qualitative (and very small quantitative) changes on the 
results on the employing of other even powers on /. 

The relatively simple form of So and ^^^ allows, in prin- 
ciple, the Hamiltonian to be evaluated analytically. How- 
ever, the number of terms involved are of the order of 
500, and even prohibitively higher when a larger number 
of excitations is introduced. Fortunately, the localized 
nature and smoothness of the excitation allows one to 



construct accurate algorithms for the evaluation of H. 
For the evaluation of the integrals on Eq. (|2]) , we used the 
composite Simpson's rule for two-dimensional integrals 
with 500 X 500 initial sampled points, and convergence 
on the 8th decimal place. The estimated maximum error 
is on the 5th decimal place (typically the 7th significant 
figure). For optimizing the free parameters of the excita- 
tion, we employed the steepest descent algorithm, with 
convergence for the energy on the 8th decimal place. 

Figure 1(b) shows the mapping of an excited state, 
S = So -\- £,re^ introduced at Xe = 7r/2, ye = in the xy 
plane, for a = 5.419 x 10-^ 7 = 0.3426, 6> = tt and A = 1. 
The parameters were obtained for A = 10~^ and /i = 2.6 
(the small value of A was chosen in order to minimize 
the string tension). We can see that the excitation fiips 
the spins only in a localized region, as expected, and the 
two adjacent vertexes no longer obey the ice rule. In the 
discrete lattice, this corresponds to the emergence of the 
magnetic monopoles. This is also the case of the continu- 
ous model - Figure 2(a) shows the density of Hamiltonian 
at the vicinity of the excitation. One can clearly distin- 
guish two localized distributions of energy, correspond- 
ing to two quasi-particles in the system. The density of 
fiux diagram in Figure 2(b) confirms that the two quasi- 
particles have opposite magnetic charge. The results are, 
obviously, dependent on the Hamiltonian coefficients A 
and /i. For /i = 0, a quick evaluation of Eq.(j2j) shows 
that the excitation is not stable - the energy monotoni- 
cally goes down to zero as we make A ^ 0. For values 
of /i higher than certain value /j^c^ the energy has a local 
minimum for some value of A ^ 0. This can be seen in 
Figure (3) - for A = 10~^ and /i = 1.2, the energy of the 
system presents only one stable minimum at the ground 
state. For /i = 2.6, however, there is a local minimum 
at A ^ 1. So, above /j^c the system presents metastable 
excited states. The inset in Figure (3) shows the stability 
of the /i = 2.6 excited state with respect to the direction 




FIG. 2: (color online) The density of energy (a) and density 
of magnetic flux (b) of the system with the introduction of 1 
excitation. 



of the excitation - the minimum energy occurs for ^ = tt, 
opposite to the direction of Sq at f^. Figure 4 shows 
that the /i x A plane is spht into two regions, one above 
jiic, where the excitation does correspond to a metastable 
state, and one below /j^c, where it does not. The model 
Hamiltonian for a continuous spin system, Eq.(|2]), pre- 
senting a non-trivial ground state that satisfies the spin 
ice rule, together with the stability of the excited state 
S^ are the key results of this paper. As a consequence 
of these features, there is the emergence of the magnetic 
monopoles as localized energy distributions with oppo- 
site magnetic fiux. In what follows, we shall see that, as 
in the discrete 2D artificial spin ice case, the energy of 
the system grows with the monopoles separation. Also, 
we will show that, due to the overlapping of the exci- 
tations, the model presents also dipoles interaction. So, 
from a classical local model, we obtain effective nonlocal 
interactions between the monopoles or the dipoles. 

The introduction of a sequence of N excitations in the 
system at subsequent positions {vi}^ spatially separating 
the monopoles, is obtained by adding to the ground state 



the function ^^^^^^(r), where ^i(r) = e 



-'y{r—ri)^ 
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fi)^ — 2){cos{0i)i + sin{Oi)j). We are already assuming 
A = 1 and 6i giving the opposite direction of the field 
So{fi). The free parameters to be optimized are 7^ and 
ai^ i = 1..N. Symmetry can be used to reduce the num- 
ber of parameters, so, for the introduction of two adjacent 
excitations, only one set of parameters is necessary. For 
a large enough in-line chain of excitations, the parame- 
ters for the excitations at the middle of the chain do not 
differs significantly from each other, so, in general, only 
four sets of parameters are necessary. As in real spin ice 
systems, this string is not unique, and, as in Refs.[l|, [9|, 
the results presented here for the system energy are for 
the shortest string connecting the two monopoles. 

The energy of the system as a function of the 
monopoles separation r is presented in Figure 5. As we 
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FIG. 3: The system energy as a function of the excitation 
amplitude A. The absolute minimum occurs at A = 0, the 
ground state. For /a > /j^c, there is a local minimum at A ^ 1. 
On the inset - The energy as a function of the excitation 
direction. The local minimum occurs when the excitation 
direction is opposite to the ground state direction. 



can see, the energy grows with r, an evidence of the pres- 
ence of the string joining the monopoles, as in the discrete 
case [9]. From N > 4 approximately, the energy grows 
quite hnearly with r. As a consequence, the monopoles 
are confined, i.e., from a specific r value it would be en- 
ergetically more favorable to disrupt the string in order 
to form two new isolated dipoles. However, at zero tem- 
perature the strings are stable - there are no fiuctuations 
that could bring the large dipole configuration to the two 
smaller dipoles one. In actual 2D spin ice systems, the 
temperature necessary to fiip a magnetic domain is pro- 
hibitively high, so there is no possibility of such a thermo- 
dynamical rupture to occu r [l3| . However, a stochastic 
external magnetic field p^6l - [l8| could provide the fiuctua- 
tions necessary to the create an effective thermal ensem- 
ble, allowing the string rupture. 

At small r, the behavior of the curve indicates that 
negative powers of r should dominate in this region, as 
also observed in ref. [9]. The dashed line in Fig. (5) rep- 
resents the least square fitting of the energy as a function 
of r to the function a-\- Pr-\-S/r. The fitting is not exact, 
but it clearly shows that there are contributions of nega- 
tive powers of r. The density of Hamiltonian for a specific 
configuration with 8 excitations is presented in Figure 6. 
We can clearly see two higher peaks, corresponding to 
the magnetic monopoles, joined by a "string" of small 
fiuctuations of the energy. 

By introducing two different chains of excitations of 
fixed length (two emergent dipoles) we can also evalu- 
ate the energy behavior as we separate the dipoles. The 
energy rapidly drops down to the noninteracting dipoles 
value {2Eo = 2 x 11.25188), since the excitations are 
characterized by a short excitation width 7 ^ 0.37. We 
evaluate the energy of a configuration of two dipoles of 
length >/27r separated by the distances r = tt and r = 27r 
in the y direction. The excitations are introduced at 
ni = (7r/2)z, ri2 = TTZ + (7r/2)j, f2i = (7r/2)i + nirj and 
r22 = TTZ + (n + l/2)7rj, with n = 1 for r = tt or n = 2 
for r = 27r. The r = tt configuration corresponds to two 
anti-parallel dipoles, and the configuration energy results 
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FIG. 4: The met ast ability diagram of the excited state. For 
the values of A and /i below the line, the excited state is 
unstable, whereas for the values above, the excited state is 
metastable. 



E = 19.54, while for r = 27r the dipoles are parahel, and 
the configuration energy is E = 22.65. As expected for a 
dipole interaction, the interaction is attractive when the 
dipoles are anti-parallel {E < 2Eq) and repulsive for the 
dipoles in the parallel configuration {E > 2Eq). How- 
ever, it is important to note that the interaction between 
dipoles depends also on the specific path employed to 
construct the strings, leading to a very rich phenomena 
to be studied. The short range dipoles interaction can 
be observed, as the configuration energy for r = Stt and 
r = 47r both agree to 2Eq at the 3rd decimal place. 

The action of an uniform external magnetic field on the 
ground state and monopoles pairs configurations can be 
easily evaluated by adding the interaction Hamiltonian 
Hh = —h. J Sdxdy to the model Hamiltonian, Eq.(|2]). 
With a not strong enough /i, in order to not change the 
system configuration, the (very simple) evaluation of Hh 
in the ground state reveals that there is no changes on 
the system energy - there is no monopoles to couple with 
the external field. With the external field aligned on the 
direction of S at the vicinity of an excitation, the sys- 
tem energy is lowered, whereas with the external field on 
the opposite direction, the energy is increased, as can be 
easily inferred from the observation of the excited state 
configuration. Fig. 1(b). So there is an effective coupling 
of the magnetic monopoles with the external magnetic 
field applied. A more interesting study to be addressed 
is how the external field can act on the present model 
in order to provide a transition of the ground state to 
another phase, and the evaluation of the influence of the 
external field on the stability of the excited solutions. 

We have presented a classical continuous spin model 
whose ground state presents 2D spin ice structure. We 
have studied the metastability of the excitations intro- 
duced in order to locally flip the spins, with the conse- 
quent emergence of magnetic monopoles. The effective 
interaction between monopoles or dipoles is analyzed. 
As a general conclusion, we obtain effective non-locality 
(among monopoles and dipoles) from a classical local 
model. Further developments of the model presented 



here are its natural extension to represent 3D spin ice 
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FIG. 5: The system energy as a function of the monopoles 
separation. 

systems and the dynamics and thermodynamics of the 
model. 
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FIG. 6: (color online) The density of energy for 8 excitation 
introduced in line. The two higher peaks correspond to the 
two monopoles. The small oscillations corresponds to the 
string joining the two monopoles. 
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